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We propose an information-theoretic statistical model to describe the universal features of those 
chaotic scattering processes characterized by a prompt and an equilibrated component. The model, 
introduced in the past in nuclear physics, incorporates the average value of the scattering matrix 
to describe the prompt processes, and satisfies the requirements of flux conservation, causality, and 
ergodicity. We show that the model successfully describes electronic transport through chaotic 
quantum dots. The predicted distribution of the conductance may show a remarkable two-peak 
structure. 
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The scattering of waves by complex systems has been 
a problem of longstanding interest in physics. 

Scattering of waves by a disordered medium has been 
studied, for instance, in optics for a long time Q. Inter- 
est in this problem has been revived, for electromagnetic 
waves and for electrons, in relation with the phenomenon 
of localization, with a host of exciting new features [§,§. 
Here, the diffusion time throughout the medium is the 
single important characteristic time. 

Examples of quantum-mechanical scattering by com- 
plex systems can also be found in nuclear and molecular 
physics. It is amazing that one can describe the scat- 
tering of a nucleon by an atomic nucleus — a complicated 
many-body problem — in terms of two distinct time scales: 
a prompt response arising from direct processes and a 
time-delayed one arising from the formation of an equili- 
brated compound nucleus. The prompt response is slowly 
varying in energy and is described by the energy aver- 
aged, or optical, scattering amplitudes; the equilibrated 
response is the difference from this energy average and is 
amenable to a statistical analysis [Q. Further examples 
of physical processes of this type, and studied in terms of 
similar notions, can be found in molecular physics, with 
interesting applications to chemistry J5j. 

Most remarkably, features similar to those appearing in 
scattering from nuclei also occur in quantum-mechanical 
scattering from simple one-particle systems |^^] . An ex- 
ample is a particle scattering from a cavity of dimensions 
larger than the wavelength, in which the classical dynam- 
ics is chaotic. One experimental realization of such sys- 
tems are the ballistic quantum dots > microstruc- 
tures in which both the phase-coherence length and the 
elastic mean free path exceed the system dimensions; the 
dot acts as a resonant cavity and the leads as electron 
waveguides. Previous work on these systems has implic- 
itly assumed the absence of any prompt response. 

Our purpose is to propose a model describing the uni- 



versal features that appear in any chaotic scattering pro- 
cess involving a prompt and an equilibrated component. 
This model was introduced in the past in the context of 
nuclear physics, using an information-theoretic approach 
pd| based on the mathematical development in Ref. [|l2| . 
Here we show that the same theoretical framework is suc- 
cessful in the description of electronic transport through 
ballistic quantum dots. We build on previous work which 
used this model in describing chaotic scattering JO and 
in simulating phase- breaking in quantum dots [|l4| . 

Semiclassical field-theoretic fllql , and random- 

matrix |ll|[l(||l7j approaches have been used to describe 
quantum transport through ballistic quantum dots. In 
Refs. the statistics of the problem was described 

by assigning to the quantum scattering matrix S an 
"equal a priori distribution" , consistent with the symme- 
try requirements. This "invariant measure" |l8| defines 
the "circular ensemble" of S matrices. The results for the 
ensemble average, variance and probability density of the 
conductance were found, in Ref. |l6| , to agree with a sta- 
tistical analysis of the numerically obtained conductance 
of a chaotic cavity connected to two waveguides, sampled 
along the energy axis. It was assumed that one could 
neglect "direct processes" caused by short trajectories 
that would give a nonvanishing energy averaged, or opti- 
cal, ^-matrix; this was enforced in the simulations with 
two stoppers which block direct transmission between the 
leads and whispering-gallery trajectories. Although the 
general problem may contain a range of relevant "time 
delays" , we present below an improvement on the above 
model in terms of two very different time scales, in a vein 
similar to the nuclear scattering problem above. We first 
summarize the information-theoretic approach |Ilf| . 

Information- Theoretic Approach — A quantum scatter- 
ing problem is described by its S matrix, which, for scat- 
tering involving two leads, each with width W and N 
transverse modes or channels, is n = 2iV-dimensional 



1 



and has the structure 
S = 



r a 

t r' 



(1) 



Here, r, t are the N x N reflection and transmission ma- 
trices for incidence from the left and r' , t' from the right. 
Current conservation requires S to be unitary, SS' = I; 
for time-reversal symmetry (as is realized in the absence 
of a magnetic field) and no spin, S is symmetric. 

Our starting point is dfj,^(S), the invariant measure 
under the symmetry operation for the universality class (3 
in question. The operation is S' = UqSVq, where Uq, Vq 
are arbitrary fixed unitary matrices in the case of unitary 
S matrices [the circular unitary ensemble (J3 — 2)], with 
the restriction Vo = Uq in the case of unitary symmet- 
ric S matrices [the circular orthogonal ensemble ((3 = 1)]. 
dfj,^\S) can be written explicitly in several different rep- 
resentations niMzi. 

The ensemble average of S, and hence the prompt com- 
ponent, vanishes when evaluated with the invariant mea- 
sure. Ensembles in which (S) is nonzero contain more 
information than the circular ensembles; they are con- 
structed by multiplying the invariant measure by a func- 
tion of S to give the differential probability 



dP, 



(0)/ 



(2) 



The information X associated with the above probability 
distribution is defined as [n9| 



J [P<S)] = / P(S){S) ln[p (S ) (S)] dn{S) 



(3) 



Far from channel thresholds, the S'-matrix is analytic 
in the upper half of the complex-energy-plane (causality) . 
We also require that the ensemble be ergodic [M , so that 
energy averages can be replaced by ensemble averages. 
These analyticity-ergodicity requirements (AE) imply the 
reproducing property 



(4) 



f((S})= / f(S)dP {s) (S) 



for a function f(S) analytic in its argument (expandable 
in a power series in S not involving S* ) . The probability 
density known as Poisson's kernel, 

Jfii (8) v - Adet(I-(S)(SW 0n+2 - W2 (5) 



I det(/-5(5t)) 



where Vp is a normalization constant, satisfies the repro- 
ducing property Eq. (Q) jl^] , and the associated informa- 
tion is less than or equal to that of any other probability 
density satisfying the AE requirements for the same (S) 
inj . Recently, this probability density for the S'-matrix 
has been derived from a statistical distribution for the 
Hamiltonian E| , explaining the coincidence noticed be- 
tween the two approaches flff]. Thus, Poisson's kernel 



describes those physical situations in which (a) the de- 
tails are irrelevant except for the average S-matrix and 
(b) the requirements of flux conservation, time-reversal 
invariance (when applicable), and AE must be met. 

For n — f, when the system is a cavity connected to 
the outside by only a one-mode lead and S — e l6 de- 
scribes reflection back into the same lead, the ensemble 
is uniquely determined by AE and a specified value of (S) 
and is given by Poisson's kernel p2j. For n > I the ad- 
ditional minimum-information criterion explained above 
is needed to determine the ensemble. 

Transport through Quantum Dots — In terms of the S 
matrix, the conductance for spinless particles is || 



G= (e 2 /h)T= (e 2 //i)Tr[tt t ]. 



(6) 



Thus, in applying Poisson's kernel Eq. ^ to electronic 
transport through quantum dots, we need the probability 
distribution w(T) of the total transmission T; i.e. 

w(T) =JS(T- Tr[i<t]) p g(S) d^(S) . (7) 

We study below the case n = 2, for which there is only 
one mode in each lead. 

We first discuss (3 — 2. Denote the elements of (S) by 



(S) = 



z y 



(8) 



where w, x, y, z are complex numbers with X = \x\ < 1, 
etc. If there is no prompt transmission, w = z = 0, we 
can perform the integrations analytically, yielding 

iu(T) = (l- X 2 f{l-Y 2 ) 2 
x{{1 - X 4 Y 4 )(\ - X 2 Y 2 ) 

- (I - T)[(X 2 + Y 2 )(l - 6X 2 Y 2 + A 4 F 4 ) 

+ 4X 2 Y 2 (1 + X 2 Y 2 )} 
+ (I - T) 2 [(I + X 2 Y 2 ){6X 2 Y 2 - X 4 - Y 4 ) 

-AX 2 Y 2 {X 2 + Y 2 )] 
+ (I - Tf(X 2 + Y 2 )(X 2 - Y 2 ) 2 } 
x{(l - X 2 Y 2 ) 2 - 2(1 - T)[(l + X 2 Y 2 )(X 2 + Y 2 ) 
- 4X 2 Y 2 } 

+ {1-T) 2 (X 2 -Y 2 ) 2 }- 5 / 2 . (9) 
For x = y the above result reduces to 



w(T) = (l-X 2 ) 

(1 - X 4 ) 2 + 2A 2 (1 + X 4 )T + 4X 4 T 2 
X [(I — X 2 ) 2 + 4A 2 T] 5 / 2 



(10) 



Calling r = 1 — X 2 , Eq. @ gives, in the limit T « I, 
the result of Ref . [Q . In the opposite case of no prompt 
component to the reflection, x = y = 0, w(T) is related 
to that of Eq. M by the replacements x — > w, y — > z, 
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T — > 1 — T. For nonzero x, y, w, z, we can express the 
result in terms of a single angular integration; it will not 
be given here because of lack of space. 

We now discuss the case /3 = 1. When (S) is diagonal 
and y = 0, one finds (see Eq. (5.13) of Ref. ]TT[] ) 



w(T) 



(1-X 2 f/ 2 .3 3 



2-y/T 



F(~,~;l;(l-T)X< 



(11) 



where F is a hypergeometric function. When (S) is di- 
agonal and x — y, we find 



,(T) 



£.( _ 

y/T \ [1 + (1 - T)X 2 - 2^1~TX cos V'] 3 / 2 



where (}^p indicates an average over -0 £ [0, 2-7r] and 



E 2 = X 2 l-T + X 2 -2VT^TXcosyj 

1-TX 2 +X 2 -2VT~^ TXcostp' 

When (S) is offdiagonal (x = y = 0) with w = z, we find 



io(T)= (1/2)(1 - Z 2 fT- 1/2 ((l + 2(2T - 1)Z 2 + Z 4 



-4VTZ(1 + Z 2 ) cos V + 4Z 2 cos 2 V)" 



•3/2 



(14) 



For several cases, these complicated distributions are 
plotted in Fig. 1 and will be discussed in connection with 
the numerical results below. 

Numerical Results — We have computed the conduc- 
tance for several stadium billiards, sketched in Fig. 1, us- 
ing the methods of Ref. p4| . w(T) was found by sampling 
in an energy window much larger than the energy corre- 
lation length but smaller than the interval over which 
the prompt response changes (so that "stationarity" , a 
condition for ergodicity M), is attained) and by using 
several slightly different structures. Typically we used 
200 energies in kW/n E [1.6, 1.8] and 10 structures found 
by changing the height or angle of the convex "bumper" 
in Fig. 1 . Thus we rely on ergodicity to compare the nu- 
merical distributions to the ensemble averages of random- 
matrix theory. In each case the optical S'-matrix was ex- 
tracted directly from the numerical data; in this sense 
the theoretical curves shown below are parameter free. 

We first consider a simple half-stadium with collinear 
leads at low magnetic field (BA/<j)Q — 2, r c ~ 55 W, 
where A is the area of the cavity, r c is the cyclotron ra- 
dius, and W is the width of the leads): w(T) is nearly 
uniform [Fig. 1(a)], and (S) is small because direct tra- 
jectories are negligible in this large structure. We thus 
obtain good agreement with the circular unitary ensem- 
ble prediction, as in previous work [ jl6| . 

In order to increase (S) we modify the situation in 
three ways: (1) introduce potential barriers at the open- 
ings of the leads into the cavity (dashed lines in struc- 
tures of Fig. 1), (2) increase the magnetic field, and (3) 



extend the leads into the cavity. The barriers (chosen 
so that the transmission of each barrier is 1/2) increase 
the direct reflection and thus skew the distribution to- 
wards small T [Fig. 1(b)]. The large magnetic field 
(BA/<j)Q — 80, r c — 1.4 W) increases one component of 
the direct transmission — the one corresponding to skip- 
ping orbits along the lower edge — and thus skews the dis- 
tribution towards large T [Fig. 1(c)]. Finally, extending 
the leads into the cavity increases the direct transmission 
in both directions and thus also skews the distribution to- 
wards large T [Fig. 1(d)]. Note the excellent agreement 
with the information-theoretic model. In panels (b)-(d) 
the curve plotted is the analytic expression of Eq. (||) 
and the corresponding one for direct transmission. 

By combining several of these modifications, different 
(S) and so different distributions can be produced. First, 
by using extended leads with barriers at their ends, one 
can cause both prompt transmission and reflection: this 
case, Fig. 1(e), is in good agreement with the predic- 
tion of the full Poisson's kernel. Finally, increasing the 
magnetic field in this structure produces a large average 
transmission and a large average reflection. The resulting 
w(T), Fig. 1(f), has a surprising two-peak structure: one 
peak near T = 1 caused by the large direct transmission 
and another near T = 1/2. Even in this unusual case, the 
prediction of the information-theoretic model is in excel- 
lent agreement with the numerical result. In these last 
two cases (e,f) the analysis was performed independently 
over four intervals of 50 energies each (since the four in- 
tervals show slightly different (S')'s ) and the four sets of 
theoretical and numerical data were then superimposed. 

Discussion — In addition to the structures shown in 
Fig. 1 we have studied cavities whose level density is 
not large enough for stationarity, and hence ergodicity, 
to hold. In this case, a sample taken across indepen- 
dent cavities for a fixed energy shows excellent agreement 
with Poisson's kernel. This suggests that the reproduc- 
ing property Eq. (|j) may be valid even in the absence of 
ergodicity; the reason for this is not understood. 

In Ref. |!| we found that increasing the magnetic flux 
through the structure beyond a few flux quanta spoiled 
the agreement with the circular ensemble; we now know 
that a nonzero (S) is generated and that the present 
model describes the data very well. The excellent agree- 
ment found here with a flux as high as 80 suggests ex- 
tending the analysis to the quantum Hall regime. 

We close by noting that the above w(T)'s should 
be experimentally accessible in structures where phase- 
breaking is small enough. Experimentally one can sam- 
ple the conductance distribution by varying the energy 
or shape of the structure with an external gate voltage 
p5[ , much as we did in collecting the numerical results. 
The barrier at the opening of the leads can be realized 
by designing a pincher gate and, of course, obtaining a 
sufficiently high magnetic field is standard. 
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FIG. 1. The distribution of the transmission coefficient for 
N = 1 in a simple half-stadium (top row) and a half-stadium 
with leads extended into the cavity (bottom row). The mag- 
nitude of the magnetic field and the presence or absence of a 
potential barrier at the entrance to the leads (marked by dot- 
ted lines in the sketches of the structures) are noted in each 
panel. Cyclotron orbits for both fields, drawn to scale, are 
shown on the left. The squares with statistical error bars are 
the numerical results; the lines are the predictions of the infor- 
mation-theoretic model, parametrized by an optical S-matrix 
extracted from the numerical data. The agreement is good in 
all cases. 



4 




Figure 1 . 



